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ABSTRACT 

We investigate the evolution of interacting binaries where the donor star is a low-mass giant 
more massive than its companion. It is usual to assume that such systems undergo common- 
envelope (CE) evolution, where the orbital energy is used to eject the donor envelope, thus 
producing a closer binary or a merger. We suggest instead that because mass transfer is super- 
Eddington even for non-compact companions, a wide range of systems avoid this type of CE 
phase. The accretion energy released in the rapid mass transfer phase unbinds a significant 
fraction of the giant's envelope, reducing the tendency to dynamical instability and merging. 
We show that our physical picture accounts for the success of empirical parametrizations of 
the outcomes of assumed CE phases. 
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1 INTRODUCTION 

Binary systems tend to shrink when they transfer mass from the 
more massive to the less massive star. The consequent contraction 
of the donor's Roche lobe leads to high transfer rates, proceeding 
on a thermal or near-dynamical timescale depending on whether 
this star is predominantly radiative or convective. There is evidently 
a tendency for such systems to enter common envelope (CE) evo- 
lution, where the donor star engulfs the system (Paczynski 1976; 
Ostriker 1976; Webbink 1984). The archetypal case arises when 
the donor star fills its Roche lobe as a giant. The strong dissipation 
caused by the accretor moving within the common envelope may 
then produce two effects: (a) the binary orbit shrinks, releasing or- 
bital binding energy, and (b) this energy release may unbind the 
envelope entirely, stopping the process before the accretor merges 
with the core of the donor. Features (a) and (b) are highly desirable 
in explaining the formation of cataclysmic variables (CVs), where 
a low-mass main sequence star orbits a white dwarf at a separa- 
tion considerably smaller than the radius of the white dwarf's giant 
progenitor, and indeed for tight binary formation in general. 

Despite this early promise, several decades of strenuous effort 
have failed to provide convincing evidence that the process works 
as originally hoped (see Taam & Sandquist 2000, and references 
therein). Technically the problem is difficult, requiring full 3D hy- 
drodynamics and a careful treatment of the dissipation processes. 
However there appear to be problems at a more basic level, which 
are manifest when we consider the most popular parametrization 
of CE (e.g. Webbink 1984). This compares the envelope binding 
energy with the change in orbital energy as 
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Here M s is the giant mass, R s its radius, M c its core mass, M c its 
envelope mass, M 2 is the accretor mass (assumed fixed during CE 
evolution) and ai, ai are the initial and final binary separations. 
Stellar structure calculations specify the dimensionless parameter 
A < 1, which is often taken as a constant, e.g. A = 0.5, and a < 1 
is a dimensionless efficiency parameter. Clearly only the combina- 
tion a\ is important. The problem which emerges is that in many 
cases the post-CE orbital separation at is too large, i.e. there was 
too little orbital energy release to unbind the giant envelope. For- 
mally this problem appears in the requirement ct\ > 1 (values 
ct\ ~ 6 are not unknown). Nelemans & Tout (2005) show that 
these problems are still worse when considering the formation of 
double white dwarf binaries, where it is possible to reconstruct the 
state of the progenitor binaries. The required values of qA are actu- 
ally negative in most cases. All these difficulties point to the need 
for an extra source of energy to unbind the giant envelope. 

A clue to this energy source comes from an observed case 
which CE evolution clearly cannot produce. The low-mass X-ray 
binary Cyg X-2 consists of a neutron star accreting from a low- 
mass donor (~ 0.5 -0.7 Mq) in a 9.8 d orbit (see King & Ritter 
1999, and references therein). The donor has radius ~ 7Rq, and 
its effective temperature implies luminosity ~ 150 Lq. Evidently 
this star was considerably more massive in the recent past and we 
are now seeing its luminous inner layers before they cool. But the 
current binary period of 9.8 d makes the binary so wide that far too 
little orbital energy was available to unbind the envelope, regardless 
of how wide the pre-CE separation ai was. However in this case the 
required extra energy source is clear. Assuming that mass transfer 
began when the donor was the more massive star and still predomi- 
nantly radiative, mass flowed towards the neutron star on a thermal 
timescale ~ 10 6 yr, giving a transfer rate > 10~ 6 Mq yr _1 . This 
is clearly highly super-Eddington. King & Ritter (1999) were able 
to suggest a plausible evolutionary sequence leading to the current 
state of Cyg X-2 by assuming that the neutron star accreted only 
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at the Eddington rate, the remaining transferred mass being ejected 
from the system with the specific angular momentum of the neutron 
star. Podsiadlowski & Rappaport (2000) similarly found a believ- 
able evolutionary history for the system starting from a somewhat 
less evolved donor. 

The crucial step in understanding Cyg X-2 is to assume that 
the accretion energy of some of the transferred matter can eject 
the remainder, without the system entering CE evolution. Although 
earlier papers considered the possibility of using accretion energy 
in this way (cf Han & Webbink 1999) they always concluded that 
the radiatively driven outflow would be dense and extended enough 
to act as a surrogate envelope for the accretor, and lead back to CE 
evolution. However King & Begelman (1999) argued that CE evo- 
lution would only occur if a sufficiently dense region of the gas flow 
on to the accretor overfilled the latter' s Roche lobe. They identified 
the boundary of this dense region as the trapping radius 
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where the accretion flow first reaches its local Eddington limit 
(Af c dd is the Eddington rate at the accretor radius i? 2 , see below) 
overfilled the accretor's Roche lobe. The trapping radius is effec- 
tively identical to the spherization radius defined for disc accretion 
by Shakura & Sunyaev (1973), who show that disc accretion at 
super-Eddington rates becomes spherical inside this radius and is 
a tenuous wind outside it, making it a reasonable estimate for the 
effective size of any dense envelope (see also section^. 

On this basis King & Begelman (1999) showed that systems 
where a neutron star or black hole accretes from a more mas- 
sive star would probably always avoid CE evolution provided that 
the donor was predominantly radiative, i.e. that mass transfer pro- 
ceeded on a thermal timescale, as in Cyg X-2. This allows values 
as high as 10 -3 Mq yr _1 , and SS433 is an example of a system 
in this state (King, Taam & Begelman 2000; Begelman, King & 
Pringle 2006). King & Begelman (1999) also suggested that such 
avoidance of CE evolution might apply even to donors which were 
largely convective, but did not attempt quantitative estimates. 

In this paper we extend the treatment of King & Begelman 
(1999) to such cases. We note that the mass transfer rates probably 
exceed the Eddington accretion rate 
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for gravitational energy release on a main sequence star (here R2 is 
the radius of the accretor, m p is the mass of a proton, c the speed 
of light and <tt the Thompson scattering cross-section). Thus it 
may be possible to avoid CE evolution even in cases where a red 
giant overflows on to a main sequence star, as required to form 
cataclysmic variables. 

Specifically we investigate here the evolution of systems in 
which mass transfer occurs from giants on to less massive compan- 
ions. We find that the trapping radius of the companion (where mass 
accretion reaches the local Eddington limit) is within its Roche 
lobe, so CE evolution may be avoided. Our treatment suggests that 
a comparison of the fractional change of mass with angular mo- 
mentum offers a more apt parametrization of the change of or- 
bital separation than the energy formalism (equation [TJ. We note 
that Nelemans & Tout (2005) have proposed such a formalism (but 
without an explicit physical reason for the mass loss) which empiri- 
cally describes the formation of sub-dwarf B star/white-dwarf and 
white-dwarf/M-dwarf binaries. We shall show that our physical 



treatment gives similar results, thus offering a possible explanation 
for its success in this case. 



2 EVOLUTIONARY METHOD 

We consider donor stars in initial evolutionary states ranging from 
the start of the Hertzsprung gap to the tip of the asymptotic giant 
branch. 

We define the mass ratio q = M1/M2 (normally > 1) with 
Mi the donor mass and M2 the accretor mass. The Roche lobe 
radius (Rl) as a function of orbital separation (a) for the donor is 
given by Eggleton (1983) 

0.49q 2/3 
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For given masses we can find the radius of the donor when it 
fills its Roche lobe at a given initial period (Pi). Using the Eggleton 
stellar evolution code (Eggleton 1971; Pols et al. 1998, and refer- 
ences therein) we then find the corresponding core mass (M c ) of 
the donor at this radius. We assume that mass transfer is so rapid 
that the core mass does not significantly increase during the mass 
transfer phase. 

The instantaneous mass transfer rate is given by tracking the 
overfilling of the donor's Roche lobe as given by Paczynski & 
Sienkiewicz (1972). Substituting their equation (4) into their equa- 
tion (A26) we find 
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where Si is 0.215; m is the fractional core-mass of the 
donor (M c /Mi); E(m) is a function tabulated by Paczynski & 
Sienkiewicz (1972); fj. is M\/{M\ +M 2 ); ARi is the difference in 
stellar radius and Roche-lobe radius (Rl) of the donor and W(/J.) 
is given by 
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We let a fraction tj) of the transferred matter be blown away (cf 
Bhattacharya & van den Heuvel 1991, who call this fraction /3). An 
expression for t[> which accounts for the proportion of matter which 
must be accreted to power the outflow is given by Han & Webbink 
(1999) (who call this quantity 1-/3) 
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where L co \d is the Eddington luminosity 
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and 0li and <f>~R2 are the potentials at the inner Lagrangian point 
and the surface of the accretor respectively and are given by 

GMi GM 2 G(Mi + M 2 ) 



a - b b 2a 3 
GMi GM 2 G(Mi + M 2 
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where b is the distance from the accretor to the inner Lagrangian 
point and is given by Warner (1976) 



- = 0.5 - 0.227 log q. 

a 



(11) 
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Figure 1. The initial periods at which a 1 Mq donor fills its Roche lobe 
and the corresponding donor core masses for various accretors: 0.2 Mq 
main-sequence stars (black circles), white dwarfs (red stars), 0.6 Mq 
main-sequence stars (green triangles) and 1 Mq main-sequence stars (blue 
squares). 



if 0.1 < q < 10 or Kopal (1959) if q > 10 

b w 2 w 3 

- = w- — — , 
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where 

If i?2 <C a equation 10 becomes for Mi > M e dd 
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This is equivalent to assuming that the transferred matter cannot be 
accreted at a rate faster than the Eddington rate (equation and 
that all transferred matter is accreted for sub-Eddington transfer 
rates. We consider below how much angular momentum this mass 
carries off. 

The trapping radius of the accretor, equation {2j, depends only 
on the mass transfer rate from the donor. During the evolution we 
compare the trapping radius with the Roche radius of the accretor 
found using equation © with the mass ratio now inverted. If the 
trapping radius is greater than the Roche lobe we assume that the 
system enters CE evolution and do not follow its evolution. In the 
rest of the paper we show only systems where the trapping radius 
never overfills the Roche lobe. 

The resultant change in orbital separation depends on the spe- 
cific angular momentum of the lost mass. If the transferred matter 
forms a disc and mass is lost with circular symmetry from this disc, 
the appropriate value is the specific angular momentum of the ac- 
cretor. As the accretor is further from the centre of mass this pro- 
cess shrinks the binary. However if the transferred matter cannot 
circularize and form a disc, i.e. its circularization radius 
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Figure 2. The initial periods at which a 1.8Mq donor fills its Roche lobe 
and the corresponding donor core masses for various accretors. Symbols are 
the same as in Fig.[T] 



(where b is the distance of the accretor from the inner Lagrange 
point Li, cf Frank, King & Raine 2002) is smaller than the trap- 
ping radius 7?tra P , the lost matter has specific angular momentum 
closer to that of L\. This follows because, unless a disc forms, the 
gas stream from the donor effectively conserves its specific angu- 
lar momentum at a value close to that of L\ (King, Whitehurst & 
Frank 1990). The donor thus loses matter with specific angular mo- 
mentum slightly less than its centre of mass, i.e. the donor slightly 
increases its specific angular momentum. Since the binary mass 
also drops, this tends to lead to expansion, as equation (20) of van 
Teeseling & King (1998) shows. This of course then makes Rd TC 
increase above i?tra P , creating a stable feedback so that the binary 
subsequently evolves with Rdrc — Rtmp- This occurs particularly 
in systems with large mass ratios, where most of the transferred 
matter is lost. These systems thus do not merge unless the trapping 
radius exceeds the Roche radius of the accretor and a CE phase 
occurs. 

To demonstrate the trapping radius is the boundary of the 
dense region we compute the density p(R) in the outflow from 
the trapping radius. This is essentially a spherical wind with 
constant velocity equal to the escape velocity v csc (Rt la , p = 
(2GM/i?trap) 1//2 from the trapping radius. Thus 
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Here we use l[3} to replace M c dd in the first factor on the rhs, giving 
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Figure 3. Final period and accretor mass versus core mass, and final versus initial period, for mass-transfer from a 1 Mq donor. Symbols as in Fig.fJ Diagonal 
lines indicate where the final and initial orbital periods are equal. 
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Figure 4. Final period and accretor mass versus core mass, and final versus initial period, for mass-transfer timescale from a 1.8 Mq donor. Symbols as in 
Fig.0 Diagonal lines indicate where the final and initial orbital periods are equal. 
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where TO2, V2 are M 2 , R2 in solar units. (Note that eqn[T8] repro- 
duces the familiar result that the optical depth through an outflow 
at the Eddington rate is of order unity at the Schwarzschild radius, 
i.e. setting M out = M ed d,R2 = iitrap = Rs = 2GM/c 2 gives 
r = pa-rRs/mj, ~ 1.) 

The density d 1 9b implies a drag force ~ irR^pv 2 , where R2, v 
are the companion's radius and orbital velocity, or equivalently a 
drag torque 



j ~ nR^apv 2 



(20) 



where a is the separation. The companion's angular momentum is 



J~(4/3)7rp 2 / 3 a 3 (a7P) 
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where p 2 is its mean density (writing its radius as fa). Using this 
in J we get the drag timescale as 



^dr; 
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In other words, the drag is slow compared with the dynamical 
timescale P provided that p < /)2. Since the companion is (at 
worst) a MS star with p 2 ~ 1 we see from J 1 9b that the drag is 
small provided a 3> f?tra P since M out S> M e dd by hypothesis. 

Once the donor fills its Roche lobe we consider the transfer of 
small portions of the envelope AM. The change in orbital separa- 
tion of the system Ao is 



Aa _ 2 Aj_ | A(Mi + M 2 ) 
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and equation d23 b becomes 

Aa 2^qAMy ipAM 
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We choose AM so that Aa is always less than 1 per cent of a. 
This formulation gives the same results as the calculations of Bhat- 
tacharya & van den Heuvel (1991) and Kalogera & Webbink (1996) 
in which the final orbital separation (af ) (after transfer of mass with 
a part ip lost with the angular momentum of the accretor) is given 
by 
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Figure 5. Figure showing the evolution of a 1 Mq donor and a 0.4 Mq 
white dwarf with an initial orbit period of 1000 days. The top panel shows 
the evolution of the radius of the donor (red line) and the donor's Roche- 
lobe (black) with donor mass. The second panel shows the masses of the 
components (blue and red for donor and accretor respectively). The third 
panel shows the mass transfer rate and the Eddington limit (red line). The 
bottom panel shows the evolution of the orbital period. 
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if ip is one. The subscript i indicates the initial value before mass 
transfer and Mu is the final donor mass assumed to be donor core 
mass upon filling its Roche lobe. The equations above are identical 
with equations (A. 5) and (A.7) of Bhattacharya & van den Heuvel 
(1991) with the substitutions ip — » /3 and appropriate notation for 
the initial and final masses (e.g. Mu — > m?). 

After each change in orbital separation we calculate the new 
donor radius, the corresponding new mass transfer rate (equation[5j 
and 7. The donor radius is found by using the radius-mass ex- 
ponents of the donor. Hjellming & Webbink (1987) calculated 
radius-mass exponents £ for polytropic donor envelopes losing 
mass rapidly. Soberman, Phinney & van den Heuvel (1997) showed 
that these can be replaced by functions of fractional core masses 
(m). We use the following prescription which is accurate to better 
than a percent 

2/ m \ 1/1-ra 
^ ~~ 3 U - m) ~ 3 U + 2m 
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Figure 6. Figure showing the evolution of a 1 Mq donor and a 0.4 Mq 
white dwarf with an initial orbit period of 100 days. The nomenclature is 
the same as Figure|5] 
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While condensed polytropic models (Hjellming 1989, 1990) gener- 
ally represent the response of red giants well, there are some cases 
where numerical models predict expansion. The resulting larger 
transfer rate and trapping radius would increase the chance of en- 
tering CE evolution. Evidently one would need to check for this in 
individual cases. 



3 RESULTS 

We consider the evolution of 1 Mq and 1.8 Mq donors transferring 
matter to four possible accretors. We choose these donors because 
they have degenerate cores on the giant branch. The accretors are 
(a) a low-mass main-sequence star with an initial mass of 0.2 Mq 
and a radius of 0.2 Rq, (b) a white dwarf of mass 0.4 Mq and a 
radius of 10000 km, (c) a main-sequence star of 0.6 Mq with an 
a radius of 0.6 Rq, and finally a solar-type star with initial mass 
1 Mq and radius 1 Rq . 

Figures[T]and|2]show the periods at which the donors fill their 
Roche lobes and their corresponding core masses. Black circles in- 
dicate the low-mass main-sequence companion, red stars the white 
dwarf, green triangles the 0.6 Mq main-sequence star and blue 
squares the Sun-like star. 

Figure [3] shows the outcomes when mass is lost from a 1 Mq 
donor. Most low-mass main-sequence stars end up with the trap- 
ping radius greater than their Roche-lobe radii and so enter CE 
evolution, producing a merged object with a portion of the donor 
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Figure 7. Figure showing the evolution of a 1 Mq donor and a 0.6 Mq 
accretor with an initial orbit period of 100 days. The nomenclature is the 
same as Figure \5\ 



envelope ejected from the system before the CE phase. For the 
white dwarf companions around 0.2 Mq is accreted during the en- 
tire evolution and the majority of this after the dynamical mass- 
transfer phase (see below). The more massive main-sequence stars 
accrete significant mass. This mass is accreted, however, after the 
dynamical mass-transfer such that the stars can cool the material 
sufficiently as they accrete it. In a number of these evolutions the 
mass-ratio is inverted during the dynamical mass-transfer phase 
which brings an abrupt end to the dynamical phase. The systems 
then evolve to longer periods (see below). 

Figure[4]shows the outcomes when mass is lost from a 1.8 Mq 
donor. Only some of the white-dwarf and main-sequence compan- 
ions survive. The white-dwarf accretors gain little mass while the 
main-sequence accretors gain significantly at the end of the dynam- 
ical mass-transfer phase (see below). We now consider the evolu- 
tion of specific systems which illustrate the general trends shown 
in these sequences. 

Figure [5] shows the evolution of a white-dwarf accretor with 
an initial orbital period of 1000 days. The donor overfills its Roche 
lobe on the asymptotic giant branch (AGB) with a core mass of 
0.56 Mq. After transferring a small portion of mass on the nuclear 
expansion timescale for the giant the donor becomes sufficiently 
oversized with respect to its Roche lobe that dynamical mass trans- 
fer occurs. During this phase 0.3 Mq of the envelope is transferred 
reaching peak mass transfer rates of just under 10 -2 Mgyr -1 . The 
white dwarf accretes at the Eddington rate during this phase and as 
it is a rapid phase (a timescale of order 100 yrs) the white dwarf 
does not gain much mass compared to its initial mass. After this 
dynamical phase however, the star has contracted back inside its 
Roche lobe and mass transfer is limited to the nuclear expansion 
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Figure 8. Figure showing the evolution of a 1.8 Mq donor and a 1.0 Mq 
accretor with an initial orbit period of 300 days. The nomenclature is the 
same as Figure \5\ 



of the envelope again. We consider this portion of the evolution to 
be conservative (all transferred mass is accreted) and the remaining 
0. 1 Mq of the envelope is accreted on to the white-dwarf. 

Figure|6]shows the evolution of a white-dwarf accretor with an 
initial orbital period of 100 days. The donor overfills its Roche lobe 
on the giant branch with a core mass of 0.35 Mq. During the dy- 
namical mass-transfer phase 0.45 Mq of the envelope is transferred 
reaching peak mass transfer rates of 10 -2 Mgyr -1 . At this mass- 
transfer rate the trapping radius (which only depends on the transfer 
rate) fills the circularization radius of the accretor. The system then 
switches to losing mass with the specific angular momentum of the 
inner Lagrangian point which tends to widen the binary (see bottom 
panel of Figure[6}. The evolution through this phase then continues 
so that .Rti-ap is the same as R c h c - If the trapping radius is larger 
than the circularization radius then the binary widens and the trans- 
fer rate drops while if is smaller the binary tightens and the transfer 
rate increases. After this dynamical phase the remaining 0. 18Mq 
of the envelope is accreted on to the white-dwarf. As the mass ratio 
is inverted the binary becomes wider. 

Once the dynamical mass-transfer phase has finished the evo- 
lution of the system will proceed much more slowly. Clearly the 
assumption that the core mass stays constant is no longer valid. The 
only effect this will have on the calculations, however, is that less 
mass is available to transfer. This will result in evolution ending 
sooner and the systems not widening as much after the dynamical 
mass-transfer stage. In the case of white-dwarf accretors the mass 
transfer has been assumed to be conservative. During the dynam- 
ical mass-transfer phase the accretion rates will be much higher 
than those required for stable nuclear burning (see Kahabka & van 
den Heuvel 1997) but afterwards may be low enough for nova be- 
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haviour. In this case the final part of the mass-transfer may be non- 
conservative and the systems would end up tighter i.e. shorter peri- 
ods. 

Figure UJ shows the evolution of a 0.6 Mq main-sequence star 
with an initial orbital period of 100 days. The donor overfills its 
Roche lobe on the giant branch with a core mass of 0.34 Mq. 
During the dynamical mass-transfer phase 0.3 Mq of the enve- 
lope is transferred reaching peak mass transfer rates of under 
10 -2 Mgyr -1 . Near the end of this dynamical phase however, the 
mass ratio has been inverted and mass transfer widens the binary 
rapidly halting the dynamical mass-transfer phase. The remaining 
0.35 Mq of the envelope is then accreted on to the main-sequence 
accretor raising its mass to over 1 Mq . The evolution of this system 
is substantially different to that of standard CE evolution where the 
accretor would be assumed not to gain in mass. 

Figure [8] shows the evolution of a binary composed of a 1.8 
and a 1.0 Mq main-sequence star with an initial orbital period of 
300 days. As in Figure [7] near the end of the dynamical phase the 
mass ratio is inverted. The system then evolves to longer periods. 

The radial response of low-mass main-sequence star to accre- 
tion is significantly different to that of a high-mass star. As Fuji- 
moto & Iben (1989) show a low-mass star with a convective en- 
velope does not significantly expand on accretion and could even 
contract. In our scenario the entropy of the transferred mass will be 
higher than the main-sequence surface but it will be subject to a far 
higher surface gravity. Regardless of the detailed radial response of 
the secondary it is clear from Figures[7jand[8]that only a relatively 
small portion of mass is accreted in the build-up to and during dy- 
namical mass-transfer and so we do not expect significantly expand 
during this process. 



3.1 Parametrizations 

Treatments of CE evolution have generally been expressed in terms 
of empirical parametrizations (Webbink 1984; Nelemans & Tout 
2005). Here we compare our results with these parametrizations. 
This serves two purposes. First, where the parametrizations are suc- 
cessful, our treatment should agree, and indeed offer an explana- 
tion for the success. Second, a simple and robust parametrization 
would give a succinct way of summarizing complex physics. As 
these parametrizations assume no accretion onto the secondary we 
make the same assumption in order to compare our model to their 
parametrizations . 

Our picture predicts a total change in angular momentum 
(A J) proportional to the total change in system mass (AM) 



AJ 



AM 



J 7 M 



(30) 



This is the form assumed empirically by Nelemans & Tout (2005). 
Our physical picture of the process allows us to specify an upper 
limit on 7 which we can compare with the results they get by ap- 
plying their parametrization to observed systems. 

Assuming that the entire envelope of the donor is lost d30t 
becomes 

AJ _ M cnv 
J ~ 1 M C + M 2 ' 
The orbital angular momentum of the system is 



(31) 



J = MiM 2 



Go 



Mi + M 2 ' 
so the left hand side of equation J30b is 



(32) 



AJ 

Ji 



M 1{ M 2{ / (Mu + M 2i )a f 



(33) 



Mu M 2i V (Mu + M 2f )oi 

If the matter is lost with the angular momentum of the accreto£] 
then taking M2 as fixed, and using equation l !28t this becomes 



AJ 

Ji 



Mu + M 2 



exp 



Mu + M 2 
Hence in this case we have 
Mu + M 2 (Mu 



( M env \ 
V M 2 J 



M 2 



7 = 



Me, 



M cl 



■ exp 



M 2 



(34) 



(35) 



,(M C + M 2 ) 

This is an upper limit as we assume that the lost matter carries 
the specific angular momentum of the accretor which will have a 
higher specific angular momentum than the inner Lagrangian point. 
Hence if there is any mass loss without a disc the specific angular 
momentum of the lost matter will be smaller and the total A J will 
be less. 

Fig.[5]shows 7 versus donor mass for various core and accretor 
masses. Nelemans & Tout (2005) show that observed close white- 
dwarf binaries are well described by the choice 7 ~ 1,5- 1.75. 
Since many of the values in Fig.|9]are within this range we conclude 
that the evolution of close white-dwarf binaries is compatible with 
the idea that all the transferred matter is ejected with the specific 
angular momentum of the accretor. 

Nelemans & Tout (2005) considered neutron-star progenitors 
in the antecedents of low-mass X-ray binaries and argued that 7 
values slightly above 1.5 allow formation. From equation d35 b we 
find that a typical system with M cnv = 6Mq, M c = 2.5 Mq, M 2 
= 1 Mq has 7 = 1.56. This is low enough not to merge, but high 
enough to form a tight binary. Again these systems are compati- 
ble with all the transferred matter being ejected with the angular 
momentum of the accretor. 

We have not considered accretion on to neutron stars but such 
accretion is likely to be super-Eddington so we would expect most 
of the mass to be ejected with the specific angular momentum of 
the neutron star. This evolution is similar to Cyg X-2. Applying 
the results of King & Ritter (1999) with Mu = 3.6 Mq , M onv = 
3 Mq , M 2 = 1.4 Mq, equation (HU gives a 7 of 1.2. 

We now turn to the comparison with the energy or aX formal- 
ism (Webbink 1984). Making the same assumptions as for equa- 
tion J35t (all mass ejected with specific angular momentum of the 
accretor) we find 



a A 



M c 



M c + M 2 / 
Mh + M 2 VMii 



exp 



2 (Mu - Me) 



Mu 

m7 



M 2 

M C M 2 (Ru/ai) 

Me) 



(36) 



2Mn (M u 

where Rn / ai is given by equation I0 and is a function of only the 
initial mass ratio (Mii/M 2 ). Fig. 1 101 shows the values of (aA) _1 
using these assumptions. The solid lines indicate the likely range 
of aX of 0.5- 1.0. Again one sees that at least some choices of a A 
produce the same outcomes as predicted by our calculations. How- 
ever, it is clear that for different accretor masses only a small range 
of donor masses are possible. This contrasts to the broad range of 
donor masses possible with the 7 prescription. Clearly 7 is not a 
constant but a free parameter which weakly constrains the possible 

1 van der Sluys, Verbunt & Pols (2006) have applied the method of Nele- 
mans & Tout (2005) to a larger number of systems and find the observed 
double white-dwarfs can be well explained by assuming that envelope ejec- 
tion occurs with the specific angular momentum of the accretor. 
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Figure 9. Values of 7 ma x for accretors ranging from 0.2 to 0.6 Mq . All matter lost from the system is assumed to carry the angular momentum of the accretor. 
Line colours indicate donor core masses of 0.2 (bottom) to 0.7 Mq (top) in increments of 0.1 Mq. Solid horizontal lines show the range for 7 suggested by 
Nelemans & Tout (2005). 



progenitors of a system compared to the aX formalism. Our calcu- 
lations may offer a physical justification for these parametrizations. 
However it is clear that neither parametrization gives a full descrip- 
tion of our results, which is perhaps not surprising. 



4 CONCLUSIONS 

We have studied binary evolution in cases which are normally 
thought to lead to common-envelope evolution. We have shown 
instead that the accretion luminosity is in many cases able to expel 
much of the super-Eddington mass transfer without catastrophic 
orbital shrinkage. Treatments of CE evolution have hitherto relied 
on empirical parametrizations (e.g. Webbink 1984; Nelemans & 
Tout 2005). In cases where these are successful in describing ob- 
served systems our theoretical treatment agrees quite well with the 
adopted parametrizations. This suggests that our ideas give a rea- 
sonable description of the prior evolution of cataclysmic variables, 
double white dwarf binaries and low-mass X-ray binaries. 



ACKNOWLEDGEMENTS 

MEB and LMD are PPARC research associates. ARK acknowl- 
edges a Royal Society-Wolfson Research Merit Award. We thank 
the referees for suggestions which considerably improved this pa- 
per. 



REFERENCES 

Begelman M. C, King A. R., Pringle J. E., 2006, MNRAS, 370, 399 
Bhattacharya D., van den Heuvel E. P. J., 1991, Phys. Rep., 203, 1 
Eggleton P. P., 1971, MNRAS, 151, 351 
Eggleton P. P., 1983, ApJ, 268, 368 

Frank J., King A., Raine D. J., 2002, Accretion Power in Astrophysics: 
Third Edition. Cambridge, UK: Cambridge University Press. 

Fujimoto M. Y., Iben I. J., 1989, ApJ, 341, 306 

Han Z., Webbink R. F., 1999, A&A, 349, L17 

Hjellming M. S., 1989, PhD thesis 

Hjellming M. S., 1990, PASP, 102, 374 

Hjellming M. S., Webbink R. E, 1987, ApJ, 318, 794 

Kahabka P., van den Heuvel E. P. J., 1997, ARA&A, 35, 69 

Kalogera V., Webbink R. E, 1996, ApJ, 458, 301 

King A. R., Begelman M. C, 1999, ApJ, 519, L169 

King A. R., Ritter H, 1999, MNRAS, 309, 253 

King A. R., Taam R. E., Begelman M. C, 2000, ApJ, 530, L25 

King A. R., Whitehurst R., Frank J., 1990, MNRAS, 244, 731 

Kopal Z., 1959, Close binary systems. The International Astrophysics Se- 
ries, London: Chapman & Hall, 1959 

Nelemans G., Tout C. A., 2005, MNRAS, 356, 753 

Ostriker J., 1976, in Eggleton P., Mitton S., Whelan J., eds, JAU Symp. 73: 
Structure and Evolution of Close Binary Systems Common Envelope 
Binaries, p. 206 

Paczyhski B., 1976, in Eggleton P., Mitton S., Whelan J., eds, JAU Symp. 
73: Structure and Evolution of Close Binary Systems Common Enve- 
lope Binaries, p. 75 

Paczynski B., Sienkiewicz R., 1972, Acta Astronomica, 22, 73 

Podsiadlowski Ph., Rappaport S., 2000, ApJ, 529, 946 



© 2006 RAS, MNRAS 000.QJ19] 



An alternative to common envelope evolution 9 





Figure 10. Values of (a A) ~~ 1 for accretors ranging from 0.2 to 0.6 Mq . All the transferred matter is assumed lost from the system carrying the specific angular 
momentum of the accretor. Line colours indicate donor core masses of 0.2 (bottom) to 0.7 Mq (top) in increments of 0. 1 Mq . Solid horizontal lines show the 
typical values assumed for aX of 0.5-1.0. 



Pols O. R., Schroder K.-R, Hurley J. R., Tout C. A., Eggleton P. P., 1998, 

MNRAS, 298, 525 
Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 

Soberman G. E., Phinney E. S., van den Heuvel E. P. J., 1997, A&A, 327, 
620 

Taam R. E., Sandquist E. L., 2000, ARA&A, 38, 1 13 

van der Sluys M. V., Verbunt E, Pols O. R., 2006, A&A in press 

lastro-ph/0610492l 
van Teeseling A., King A. R., 1998, A&A, 338, 957 

Warner B., 1976, in Eggleton P., Mitton S., Whelan J., eds, IAU Symp. 

73: Structure and Evolution of Close Binary Systems Observations of 

Dwarf Novae, p. 85 
Webbink R. F., 1984, ApJ, 277, 355 



This paper has been typeset from a TpX/ L5TgX file prepared by the author. 



© 2006 RAS, MNRAS OOO.QJH] 



